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A theory of light diffraction from planar quasicrystalline lattice with resonant scatterers is pre- 
sented. Rich structure, absent in the periodic case, is found in specular reflection spectra, and 
interpreted as a specific kind of Wood anomalies, characteristic for quasicrystals. The theory is 
applied to semiconductor quantum dots arranged in Penrose tiling. 
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I. INTRODUCTION 

Discovery of quasicrystals initiated new fields of re- 
search in solid-state photonicsiii^ These deterministic ob- 
jects allow Bragg diffraction of light, like conventional 
photonic crystals, but are not restricted by the require- 
ment of periodicity, and thus can be easier tailored to the 
desired optical properties. Such an extra degree of free- 
dom is especially important for the control of light-matter 
interaction in resonant photonic structures,^ where the 
constituent materials possess resonant excitations, like 
excitons or plasmons. For example, one-dimensional po- 
laritonic Fibonacci quasicrystal based on quantum-well 
excitons has been realized in Ref. 4, while the two- 
dimensional (2D) plasmonic deterministic aperiodic ar- 
rays of metallic nanoparticles have been fabricated in 
Ref. i. 

On the other hand, it is well known that interfer- 
ence between different processes of scattering on arbi- 
trary grating can lead to the so-called Wood anomalies 
in optical spectra.—^ Indeed, the incident plane wave can 
undergo either specular reflection or diffraction. The in- 
terference of these processes may result in intricate opti- 
cal spectra. To the best of our knowledge, no systematic 
study of Wood anomalies in quasicrystalline gratings has 
been performed yet, although their high importance for 
light transmission through aperiodic arrays of holes in 
metallic thin films has been mentioned in Ref. H|. Here 
we present a theory of light diffraction from the 2D res- 
onant photonic quasicrystals made of quantum dots and 
show, that lattice Wood anomalies of novel type, com- 
pletely absent in the periodic case, are manifested in their 
optical spectra. 

The rest of the paper is organized as follows. In Sec.lIIl 
we formulate the problem and outline the calculation 
technique. Sec. Mil presents approximate analytical re- 
sults for reflection coefficient. Results of calculation are 
discussed in Sec.[lV] Sec. [V]is reserved for conclusions. 




(c) 



FIG. 1: (Color online), (a) Canonic Penrose tiling, (b) Cal- 
culated diffraction image of this tiling. The diameter of each 
spot, located at the point corresponding to the Bragg diffrac- 
tion vector G, is proportional to the absolute value of the 
structure factor \fa\- Only the spots with |/g| > 0.15 are 
shown. Red circle indicates the value of the diffraction vector 
G* « (2n/a r ) x 1.05, where a r is the side of the rhombus in 
Penrose tiling. Filled and empty circles correspond to oppo- 
site stars ±G*e n , where n = 0. . .4. (c) Schematic illustration 
of the coupling between five wave vectors of the star G*e„ due 
to Bragg diffraction. 



incident wave propagates along z axis. The Penrose tiling 
shown in Fig. QJa) has fivefold rotational symmetry and 
can be denned as follows. First we introduce the basis of 
five vectors 



6n — 



cos 



2im 



, sin 



ra = 0...4. 



Then five sets of equidistant parallel lines r nj are de- 
fined, each set normal to the corresponding vector e n : 
= j + (2/5), where index j accepts all integer 
values. Finally, each cell in the obtained grid, bounded 
by the lines r 0j0 , r ,j +i, ■ . . r 4j4 , r 4j4+ i, is mapped to 
the point r = a r ^2 n=0 j n e n , belonging to the Penrose 
lattice with the rhombus side equal to a r . Such ap- 
proach is termed as dual multigrid technique^ It can be 
used to generate quasiperiodic lattices with arbitrary de- 
gree of rotational symmetry. Other equivalent definitions 
of the Penrose tiling are based on the cut-and-project 
technique^ The lattice structure factor 



II. PROBLEM DEFINITION AND METHOD OF 
CALCULATION 



1 N 

f^ = ^NT, c2iqri = E fe^-e hl ., 4 (1) 



The structure under consideration consists of quantum 
dots, arranged in the canonical Penrose tiling^ in the 
plane xy and embedded in the dielectric matrix. The 



is shown in Fig. [Tfb) , it consists of Bragg peaks at the 
2D diffraction vectors Gh x ...h4 — G* J2t=i h n e n , where 
G* = 47TT 2 /(5a r ) and r = (\/E+l)/2 is the golden mean. 
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Each non-zero diffraction vector G belongs either to the 
five-vector star Ge n or to the opposite five-vector star 
— Ge„, where n = . . .4. Since /_g = /gi the absolute 
value of the structure factor is the same for both stars 
and thus has the tenfold rotational symmetry. 
Electric field E satisfies wave equation 



by 



rot rot E(r) 



D(r), 



(2) 



where the displacement vector D(r) = SbE(r) + 
AitP cxc {r) includes nonresonant contribution, character- 
ized by background dielectric constant £b, and excitonic 
polarization P cxc . The material relation between the ex- 
citonic polarization and the electric field reads^ 

47rP cxc (r) = T(w) ^ $(r - a) / dV$(r' - a)E(r'), 

a J 

where the resonant factor T(u) is given by 



T( W ) 
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Equation ([3]) contains summation over all dots, centered 
at the points a and characterized by excitonic envelope 
functions <fr(r — a). Other excitonic parameters in Eq. (j4]) 
are as follows: longitudinal-transverse splitting wlt and 
Bohr radius in the corresponding bulk semiconductor as, 
resonance frequency ujq and phenomenological nonradia- 
tive damping T. In the following calculations the ex- 
citonic envelope function is taken in the Gaussian form 
$(r) =f x (ttR)~ 3 / 2 exp(— r 2 /R 2 ), where R is the char- 
acteristic radius of quantum dot. For a ls-exciton, quan- 
tized as a whole, one hasi 2 - $o = 2~ 3 / 2 /-\/7ra|. The 
results weakly depend on the specific choice of the en- 
velope function. For instance, similar equations can be 
applied to the cluster of small metallic spheres with the 
radius R <C 2irc/ui. In this case the functions $(r) are 
constant for r < R and zero for r > R and the factor 
T(u>) should be replaced by the resonant susceptibility of 
the metallic sphere near given plasmon resonance. 

Our calculation approach generalizes the methods used 
in Refs. [l3lfl3 . Electric field dependence on the co- 
ordinates x and y is described in the basis of plane 
waves. Different plane waves are coupled due to the 
Bragg diffraction. Coupling strength is determined by 
the structure factor coefficients. We keep in the plane 
wave expansions only 61 diffraction vectors G with the 
largest values of fc, shown in Fig. [TJb). Substituting 
Eq. ([3]) into Eq. (jSJ and applying Fourier transformation 
Ek = J d 3 r exp(—ikr)E(r) we obtain closed equation for 
the electric field: 



e ik ' a $ k/ E k , + El 0) 



fco r d?k' 

k k 2 -q 2 ^ J (2tt) 3 

Here q = Wy^Fb/c, $fe = $o exp(-i? 2 fc 2 /4) is the exci- 
tonic envelope in fc-space and the matrix Uk is defined 



[Uk\uf3 = d a p — 



The inhomogeneous term E^ in Eq. ([5]) describes the 
incident wave. We will distinguish in-plane and perpen- 
dicular components of all vectors and use the notation 
Q = (Q\\,Qz), Q\\ = {Qx,Q y ) ■ Introducing the struc- 
ture factor 



yy(fe-fe')c 



S 



G 



6(k» - fcf, - G)f G , (6) 



where S ~ 0.81a 2 is the mean area per lattice site in the 
Penrose tilingfii we obtain 



dk' 
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From now we restrict consideration to the case of normal 
incidence, E^(r) = e exp(ig,z). Multiplying Eq. (0 by 
&G,k z and integrating over k z we obtain a set of linear 
equations 



A G = x l gUg^2Jg-g'^g' + S G ,oe (8) 

G' 

for in-plane vectors 



Ag = 



/ dk s ®G,k z Ee,k, 



G J dk z &i*$ G ,K 

Other quantities in © are the dimensionless susceptibil- 
ity 



X 



cjo — w — ir 

and complex coefficients 
. q 



To p 7r a| WLT (<I>o) 2 e-(^ 2 / 2 
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Lg — i— erfc ( -^-y=- ) , <Zg = \fq 1 ~~-G 1 ~ . (9) 
qc V V2 / 



\ V2 J 

The complementary error function in ([9]) is defined as 

2 f x 

erfc(x) = 1 = / exp(-< 2 )di. 

V 7 *" Jo 

We note, that since erfc(a;) vanishes with asymptotics 
exp(— x 2 )/(y / 7rx) at x — > oo , both quantities Lq and 
Ac quickly decay at G > X/R assuring the convergence 
of the sum over G' in Eqs. (JSJ). After the vectors Ag are 
found from system (|8|), electric field is given by the in- 
verse Fourier transformation of Eq. (J7J . At the large dis- 
tances from the plane with quantum dots ( \z\ 3> R,l/q) 
we obtain 



E(r) 



G 



e- iqGZ E£\ 0-)--oo), 

e iqaz E W } ( z ^ + oo), 
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where p = (x, y) and the amplitudes of the reflected and 

J G allLt ' are 



transmitted waves, E, 



E'a> = U G , qG S G , E G > = e 5 G ,o + U G ^ qa S G , (10) 
S G = i x — YJg-g'Ag, . (11) 



G' 



Thus, the specular reflection coefficient is given by 
R(to) — \Sq\ 2 . Due to the fivefold rotational symmetry 
of the Penrose tiling the vector So is parallel to eo, and 
its magnitude is independent of the orientation of eo . For 
zero exciton nonradiativc decay (r = 0) the energy flux 
conservation condition along z direction reads 



|S | 2 + |l + So| 2 +2 Y, 

0<G<q 



qc_ 

<i 



\GS, 



G 



'G 



1 



(12) 

where three terms in the left hand side correspond to 
specularly reflected, transmitted and diffracted waves, re- 
spectively. Equation (|12[) can be rigorously derived from 
Eqs. (|5))- (ITTj) taking into account that ImLg = q/q G for 
G < q. 

Reflection coefficient r(ui) = So(u>) has a simple ana- 
lytic form only if the in-plane diffraction is totally ne- 
glected, i.e., only one vector G — is included in the 
plane wave expansions Eqs. ©. The result reads 



r(w) 



ir n 



u) - u> - i(r + r ) 

iq(uj )R 



(13) 



Gjq = ojq + Tq Im erfc 



V2 



Thus, the quantity To can be interpreted as exciton ra- 
diative decay (evaluated neglecting diffraction), while 
luo is the exciton resonance frequency renormalizcd by 
the interaction with light. Equation (|13l) is similar to 
the reflection coefficient from the quantum-well exciton 
resonance^ It is valid only if inter-dot distances are 
small compared to the light wavelength, a r <C 2n/q. 
If o, r > 27r/<7, the waves with non-zero in-plane wave 
vectors must be included into theoretical consideration. 
Generally it can be done only numerically. Analytical 
expression for the reflection coefficient, obtained taking 
into account the diffraction vector G = and the star 
of given diffraction vector G, is presented in the next 
Section. We note, that although the experimental real- 
ization of the coupling between the spatially separated 
quantum dots via electromagnetic field is a challenging 
task, a substantial progress has been recently achieved 
for dots in the microcavitiesi 15 ' 16 



III. REFLECTION COEFFICIENT IN THE 
TWO-STAR APPROXIMATION 

In this Section we consider general quasicrystalline 
tiling with A-fold rotational symmetry. Specular reflec- 
tion coefficient of the normally incident light is calculated 



analytically. We take into account N + 1 diffraction vec- 
tors, belonging to the two stars, namely, the trivial star 
of vector G = and the star of the given vector G* , 
including N diffraction vectors 



G h 



G* (cos n<p, sin rip) 



where n — 0...N — 1 and ip = 2tt/N. The coupling 
between the plane waves corresponding to G — and 
G = G n is described by the structure factor coefficient 
/g* ■ It is also essential to consider the coupling between 
the wave vectors within the star of vector G* . This cou- 
pling is characterized by structure factor coefficients 



fm — /G m+n -G„ 



(14) 



and shown for a Penrose tiling (where N = 5) in Fig.[IJc). 
Solid and dashed lines correspond to two possible values 
of coupling coefficients f m . Under the normal incidence 
of the wave eoexp(iqz) all the vectors A lie in the xy 
plane. The solutions of Eqs. (jHJ) can be sought in the 
form 



and 



A G „,i 
Ag„,„ 



Co 



C 2 



A = A e , 



&0,x 

eo, v 

cos(2rup)eo, x + sin(2ny)eo,y 
— cos(2nip)eo, y + sin(2ntp)eo >x 



(15) 



(16) 



The vector Ao in Eq. (|15p is parallel to eo and its mag- 
nitude is independent of polarization. The structure of 
Eq. (|16p is more complex. Let us examine it for the 
C$ v point symmetry group^ of the Penrose tiling. Both 
terms in Eq. (|16|) are transformed by symmetry opera- 
tions like vectors, and belong to the irreducible represen- 
tation Ei. They stem from the direct product D x E\, 
where D = A\ + E\ + E% is reducible representation de- 
scribing the transformation of functions 5{k\\ — G n ) and 
the representation E\ describes the transformation of the 
polarization vector components eo, x and eo, v - First and 
second terms in Eq. (|16p originate from the invariant A\ 
and irreducible representation Ei contained in D, respec- 
tively. 

The set of A G in ([S]) is characterized only by three 
unknown coefficients Ao, Co, C*2, and Eqs. ([8]) can be 
solved straightforwardly. The results are given below. 
The ratio between the two coefficients in Eq. (jT5|) equals 
to 



Co 
C 2 

where r\ = 1 — G 2 /q 2 , 



1 + 77 - 2r](A - B) 

7/ — 1 



A = X L G 
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It is convenient to introduce coefficient Q defined as 



fa- 



N-l 

E 

n=0 



A n EE QAq 



and given by 
„ N X L G ,\f G » 



I + r) - 2r)(A - B) 
1 - (r; + I) A + r,{A 2 - B 2 ) 



(17) 



In general Eq. (fT7| has two resonant frequencies corre- 
sponding to the mixed longitudinal and transverse modes 
excited within the star G n . At Bragg condition q = G* 
the coefficient rj vanishes, and only one mode is active: 

n _ N X L G *\f G *\ 2 



2(1 -A) 



Near the Bragg resonance q(ujo) = G the function L G * (ui) 
can be approximated by the following expression 



L G * 



U) 



2(w - w) 



(18) 



Neglecting nonradiative damping T we obtain the the 
eigenfrequency from the condition A = 1, result reads 



oj - r x 



UJq of 

~2T^ 



1/3 



JV-1 

E 

m— 



f m cos 2 mip. (19) 



Finally, the reflection coefficient is given by 

ir [i + e(w)] 



ir-ir L [i + g(u>)}-. 



w 



(20) 



To test Eq. [[20]) we consider the periodic quadratic lat- 
tice analyzed in Ref. [l3|. In this case f m = f G * — 1 and 
A = B = (N/2) X L G * , where N equals to 4 or 8 depend- 
ing on the value of G* . Thus, we get 



Q = 



N X L G [l - (G*) 2 /2q 2 } 
l-N X L G *[l-{G*y/2q*\ 



and 



ir 



u - ir {L a + nl g , [i - (G*)W]} - ir 



- ui 

(21) 

in agreement with Ref. Il3l . This expression has a pole 
at the frequency determined by the sum of contribu- 
tions of the stars, corresponding to the vectors G = 
and G — G* . Note that one should take into account 
the singular dependence of L G » on u, given by Eq. (|18|l 
for Bragg structure. This square root singularity gener- 
ally leads to the standard Wood anomalies in periodic 
lattices^ It turns out that for the relatively weak quan- 
tum dot exciton resonances (yTo/wo < 1), the features 
due to this singularity are not resolved in optical spec- 
tra of the periodic structure. However, in the quasicrys- 
talline system the expressions (|17p. (|20p have complex 
structure, determined by interplay of the resonances at 
the frequencies, given by Eqs. (fT5)l and (Q3 
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FIG. 2: (Color online). Specular reflection coefficient from a 
Penrose quasicrystal tuned to the Bragg condition uic = wo. 
Calculation was performed for the following set of parame- 
ters: r /w = 1CT 3 , r/r = O.l, R/a r = 0.2. Solid curve 
presents numerical calculation including all diffraction vec- 
tors from Fig.[]Jb), while dotted curve is the analytical result 
(|20p . obtained including 11 vectors. The vertical line shows 
the position cj* of the resonance of the star of vector G* , de- 
fined in Eq. (|19[) . Dashed curve corresponds to the reflection 
coefficient from a periodic square lattice, see Eq. (|21[) . Other 
parameters are indicated in text. The inset schematically il- 
lustrates the geometry of the problem. 



IV. RESULTS AND DISCUSSION 

The results of calculation for a periodic and Penrose 
lattices are presented in Fig. [2j Dashed curve shows re- 
flection coefficient for a periodic square lattice tuned to 
the Bragg resonance at the diffraction vector \G\ = G* = 
47rr 2 /(5a r ): 



uj = u G * = cG* /rib 



(22) 



This calculation was performed analytically in the two- 
star approximation Eq. (|2ip with N = 4 vectors in the 
star. In agreement with Ref. [l3l only one Lorentzian 
peak is present in reflection spectrum. The red shift of 
the peak frequency from exciton resonance luq is due to 
radiative corrections. The absolute value of the shift is 
large compared to To, because the structure is tuned to 
the Bragg condition (f2"2"|) . No Wood anomalies are re- 
solved for a periodic structure since Tq <C too- 

Dotted curve in Fig. [2] depicts analytical result (|2"0")) 
for a Penrose tiling. The spectrum was calculated taking 
into account the diffraction vector G = and ten vectors 
±G*e„, belonging to the two opposite stars. These stars 
are not coupled by diffraction and interact with the wave 
G = independently. Absolute values of the structure 
factor coefficients are the same for both stars. Thus, the 
two-star approximation presented in Sec. lIIII can be easily 
extended to include three stars G = and G = ±G*e n 
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FIG. 3: (Color online) Specular reflection coefficient from 
Penrose quasicrystals with different lattice constants a r , cal- 
culated for lug* /wo = 0.985, 1, 1.05, 1.2. For better presen- 
tation each following curve is shifted upwards by 0.4 from 
the preceding one. Horizontal lines mark corresponding zero 
levels. Other parameters are the same as for Fig. [21 



in Fig. [2j This splitting is on the order of Tq^o/To) 1 / 3 
and can considerably exceed Tq. Such spectral structure 
may be also observed in the slightly distorted 2D periodic 
resonant structure^ and in the 2D resonant photonic 
crystals with compound elementary supercellJ^ To ob- 
tain the spectra with the shape similar to those in Fig. [5] 
it suffices to tune the structure to the Bragg condition 
with the structure factor coefficient large but less than 
unity. 

Fig. [3] illustrates effect of detuning from the reso- 
nant Bragg condition (|2"2")l . One can see that the split- 
ting increases with detuning for loq* < ujq. The low- 
frequency peak becomes sharper and its position for the 
large detuning, |w — u)g*\ 3> Fo^o/Fo) 1 / 3 , is close to 
loq* . For ojg* > wo the value of the splitting decreases. 
With the further increase of detuning the peaks merge 
into one peak with structured dips, similarly as it oc- 
curs in the one-dimensional Fibonacci multiple quantum 
wellsJ^ A qualitative difference between spectral shapes 
for ujq* < wo and u>g* > Wo, revealed in Fig. [3] can be un- 
derstood taking into account that the diffraction channel 
with G = G* is closed for u)q* > w , i.e., the correspond- 
ing wave exp(i(7(3* \z\) becomes evanescent. 



V. CONCLUSIONS 



by making the replacement Q — > 2Q in Eq. (|2"TJ)) . The 
values of structure factors used for analytical calculation 
are \f G *\ = 0.38, f ±1 = 0.47 [solid lines in Fig. [He)], 
/ ±2 = 0.76 [dashed lines in Fig. (He)]. 

Allowance for non-zero diffraction vectors leads to the 
splitting of the single resonance ([T3"|) into two peaks. The 
resulting spectrum remarkably differs from that in the 
periodic case, cf. dashed and dotted curves in Fig. [21 
and is well described by the chosen approximation. In- 
deed, since the structure is tuned to the specific Bragg 
resonance ([22]) . the spectrum does not change consider- 
ably when extra diffraction vectors with G ^ 0, G / 
G* are taken into account (solid curve). Thus, Fig. [5] 
demonstrates, that the complex structure of the reflec- 
tivity spectrum due to the multiple grating resonances is 
the characteristic property of the quasicrystalline lattice. 
The magnitude of the splitting is approximately given by 
the difference between ujq and the position of the reso- 
nance w*, given by Eq. (fT9| and denoted by vertical line 



To summarize, we have developed a theory of light 
diffraction on the 2D quasicrystalline planar array of 
quantum dots. An analytic expression for the reflection 
coefficient has been derived. While, for the periodic lat- 
tice, the specular reflection spectrum has a single peak 
near the exciton resonance frequency, for the quasicrys- 
talline lattice the spectrum consists of two peaks. The 
more complicated structure of the spectrum is related 
to the interplay between specular reflection and in-plane 
light diffraction, resulting in Wood lattice anomalies, spe- 
cific for the resonant quasicrystalline structures. 
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